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Size-dependent rheology of type-l collagen networks 

Richard C. Arevalo, Jeffrey S. Urbach and Daniel L. Blair 
Department of Physics, Georgetown University, Washington, DC, 20057, USA 

ABSTRACT We investigate the system size dependent rheological response of branched type I collagen gels. When sub- 
jected to a shear strain, the highly interconnected mesh dynamically reorients, resulting in overall stiffening of the network. 
When a continuous shear strain is applied to a collagen network, we observe that the local apparent modulus, in the strain- 
stiffening regime, is strongly dependent on the gel thickness. In addition, we demonstrate that the overall network failure is 
determined by the ratio of the gel thickness to the mesh size. These findings have broad implications for cell-matrix interac- 
tions, the interpretation of rheological tissue data, and the engineering of biomimetic scaffolds. 
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The increasing variety and availability of purified extra- 
and intra-cellular matrix proteins provides an unprecedented 
opportunity to quantify the mechanical properties of the cel- 
lular environment. Reconstituted biopolymer fiber networks 
can exhibit nonlinear rheology that is dramatically different 
than synthetic polymer gels [1, 2, 3, 4, 5]. Tying the micro- 
scopic physical behavior of the network constituents to the 
meso- and macro- scale rheology is critical for determining 
the mechanosensory mechanisms that underlie such diverse 
processes as cell motility, cancer metastasis and tumor pro- 
liferation. In vivo, cells sense, and respond to forces gener- 
ated within or transmitted through the extracellular matrix 
(ECM). In turn, the regulation of sensory cues can be highly 
susceptible to changes in the ECM stiffness on length scales 
comparable to the cell size. However, in spite of the increas- 
ing number of experimental and theoretical investigations of 
the nonlinear rheological behavior of semi-flexible and stiff 
polymer networks [6, 7, 8, 9], the effects of system size on 
strain-stiffening remains essentially unexplored. Most theo- 
retical models focus on the global response of biopolymer 
networks using three distinct microscopic mechanisms; non- 
affine deformations of crosslinked rigid rods [9], entropic 
penalties of fiber stretching in entangled networks [3], or 
changes to the network geometry [10]. There is little ex- 
perimental data to directly verify these microscopic mech- 
anisms, in part because of the difficulty in extracting the 
relevant quantities experimentally [11, 12]. In addition, the 
nonlinear rheology of stiff networks can only be investigated 
with bulk techniques, so the behavior at biologically relevant 
microscopic and mesoscopic scales is largely unknown. 

Collagen is the most abundant ECM protein, and a com- 
mon component of in vitro cell cultures and bioengineered 



scaffolds. Under appropriate conditions, ry/?^ / collagen self- 
assembles in vitro to form percolated networks (gels) [13] 
with mesh sizes ^ that depend on the concentration and poly- 
merization conditions (Fig. lA-C) [14, 12]. These branched 
networks substantially strain-stiffen over a broad range of 
concentrations [15, 16, 17]; one interpretation of the biome- 
chanical function associated with nonlinear stiffening is a 
passive protection for soft tissues when shear deformations 
become large. Collagen networks are not only excellent 
ECM model systems, but provide a platform for investigat- 
ing cellular motility within three-dimensional microenviron- 
ments [18], and provide bio-compatible scaffolds for tissue 
growth and organ regeneration [19, 20, 21, 22]. 

In this work, we investigate strain- stiffening of collagen 
networks under steady shear and observe that the nonlin- 
ear rheological response is strongly dependent on the ma- 
terial thickness. Moreover, the apparent moduli near yield 
decreases dramatically in thin samples in a process that is 
controlled by the ratio of the sample size to the network mesh 
size. The system size dependence is not accounted for in cur- 
rent models of nonlinear strain- stiffening in biopolymer net- 
works. 

We apply continuous shear strains to collagen networks 
using an Anton-Paar MCR-301 bulk rheometer with a 25 
mm diameter parallel-plate geometry within a temperature 
and humidity controlled environment. All continuous shear 
strain experiments are performed at a strain rate of 7 = 1.0% 
s~^. Type Irat tail collagen (BD Bioscience, San Mateo, CA, 
3.27 mg/mL) is polymerized at 23^C for 45 minutes in lOX 
PBS at pH 7.0 with ionic strengths / = {0.044, 0.087, 0.13} 
for the concentrations c = {1.0,2.0,3.0} mg/ml, respec- 
tively. 
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FIGURE 1 Confocal reflectance images of branched type I rat 
tail collagen fiber networks corresponding to the concentra- 
tions: (A) 1 mg/ml, (B) 2 mg/ml, (C) 3 mg/ml. Each image is a 
10 /im thick maximum-projection; scale bar = 25 /im. 



The plate tool provides precise control of the sample thick- 
ness through a change of the rheometer gap h. At low strains, 
the measured stress a is proportional to the applied strain 7, 
for all values of h. Intermediate strains reveal that the net- 
work undergoes a substantial nonlinear increase of a, indica- 
tive of strain- stiffening [3]. Specifically, for h > 150/im, 
(7(7) is /i-independent. However, for smaller gaps, we ob- 
serve that the maximum of (7(7) moves to larger values of 7 
(Fig. 2). The magnitude of a at the peak remains relatively 
constant. The yield strain in thin samples is more than twice 
that of thick samples, and the apparent modulus of thin sam- 
ples can be up to 3 times smaller than the bulk value. To 
ensure that the changes in h are not introducing variability 
in the rheological measurements, particularly at small h, we 
measure the linear viscoelastic modulus Go using oscilla- 
tory rheology at a frequency of = 27r rads s~^ at a strain 
of 7 = 1.0%. For each concentration, we observe that Go is 
constant over the entire range of h (Fig. 3A). Higher strains 
lead to irreversible deformations and yielding, as evidenced 
by a reduction in the overall modulus. In addition to this dis- 
tinctive rheological signature, we observe a dramatic change 
in the rheological response by varying the sample thickness 
from h = 50 — 300/im, at a fixed collagen concentration 
(c= 1.0 mg/ml). 

To quantify the /i-dependence of the nonlinear rheology, 
we perform a least squares polynomial fit to the peaks of 
the rheology curves in Figure 2A. We use the peak posi- 
tions to define the yield values for the stress ay and strain 
jy. We determine jy as a function of h for three different 
collagen concentrations. For each concentration, 7^ exhibits 
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FIGURE 2 Stress-strain rheology curves for 1 mg/ml collagen 
gels at different gaps h (strain rate = 1% s~^); gel networks are 
polymerized at 23°C and pH 7.0. All lines connecting points are 
guides for the eye. 
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FIGURE 3 For all graphs, (o) 1 mg/ml, (□) 2 mg/ml, (A) 3 mg/ml. 
(A) Linear elastic modulus Go from oscillatory rheology versus 
gap height h. (B) Yield stress extracted from the stress-strain 
curves ay versus h shows strong ^-dependence and limited h- 
dependence for all concentrations. (C) Yield strain 7^ versus h 
for each collagen concentration extracted from the peaks in the 
stress-strain curves (Fig. 2). Error bars represent standard error 
from four trials, lines are guides to the eye. 

a clear increase for small h and approaches a constant value 
for large h (Fig. 3C). Interestingly, for the entire range of h 
and concentration the values of ay show weak or no gap de- 
pendence (Fig. 3B). We speculate that ay is determined by 
the strength of fiber- fiber or fiber-boundary junctions, and 
depends only on the total applied stress, independent of h. 

To elucidate the role of the microscopic length scale ^ in 
the observed size-dependent rheology, we replot the data in 
Figure 3C with h rescaled by ^. To quantify the mesh size ^, 
we analyze the spacing between fibers extracted along hor- 
izontal and vertical lines within z-resolved planar confocal 
images [14]. We fit the distribution of the fiber spacing to an 
exponential decay for each concentration, c = {1.0, 2.0, 3.0} 
mg/ml, to obtain the characteristic mesh sizes ^ = {14.7, 
10.4, 9.4} /im (Fig. 4A). To account for the concentration 
dependent modulus, we also rescale jy by the ^-dependent 
plateau values 7^, determined by fitting jyih) to an expo- 
nential decay with an offset. The data from the three concen- 
trations collapse onto a single curve (Fig. 4B). 

We conclude that the nonlinear rheology of collagen is de- 
termined by the interplay between the macroscopic material 
thickness and the mesh size. We do not claim that the mesh 
size determines the network rheology; this is most likely de- 
termined by some complex interplay of the network structure 
and topology. In addition, we note that the range of gaps in- 
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FIGURE 4 (A) Mesh size distributions for (o) 1 mg/ml, (□) 2 
mg/ml, (A) 3 mg/ml. Exponential fits, given by the dashed lines, 
provide characteristic mesh sizes: (o) 14.7 /im, (□) 10.4 /im, (A) 
9.4 /im. These curves have been shifted along the ordinate for 
clarity. (B) Rescaling by the extrapolated plateau yield strain 
7p and h by ^ reduces Figure 3C to a universal yielding curve 
for all concentrations; (o) 1 mg/ml, (□) 2 mg/ml, (A) 3 mg/ml. 

vestigated directly corresponds to the radial gap separation in 
many cone-plate tool geometries, and is equivalent to length 
scales of many important biological processes. 

The observed size dependence of the nonlinear rheology 
indicates that the material response is non-uniform on scales 
significantly larger than the mesh size, and is therefore in- 
consistent with models that attribute strain stiffening to the 
affine deformation of networks of semi-flexible filaments. 
Collagen fibers are relatively stiff, and therefore more likely 
described by athermal models that attribute strain stiffening 
to non-affine filament stretching [23]. The non-affine length 
scale can be significantly larger than the mesh size [24], 
but to our knowledge the implications for system size de- 
pendence has not been investigated. Alternatively, hetero- 
geneous localization of strain within the network that is ei- 
ther uniformly distributed through the gel or localized near 
the boundaries could produce size dependent effects. Re- 
flectance images do not show any measurable variation in 
collagen structure near the boundary, but more subtle bound- 
ary effects may also play a role. Finally, the long persistence 
length of the collagen fibers may directly contribute to a fi- 
nite size effect. Whatever the physical mechanism, the size- 
dependent rheology extends to important length scales for a 
range of fundamental studies and applications. 
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